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Abstract 



Recently, a new framework for describing the multiverse has been proposed which is based 
^1^' on the principles of quantum mechanics. The framework allows for well-defined predictions, 

both regarding global properties of the universe and outcomes of particular experiments, ac- 
cording to a single probability formula. This provides complete unification of the eternally 
^^..^ inflating multiverse and many worlds in quantum mechanics. In this paper we elucidate how 

^ ■ cosmological parameters can be calculated in this framework, and study the probability dis- 

^O , tribution for the value of the cosmological constant. We consider both positive and negative 

1^ I values, and find that the observed value is consistent with the calculated distribution at an 

^T) • order of magnitude level. In particular, in contrast to the case of earlier measure proposals, our 

r — I . framework prefers a positive cosmological constant over a negative one. These results depend 

^D I only moderately on how we model galaxy formation and life evolution therein. 



X 



1 Introduction 

An explanation of a small but nonzero cosmological constant is one of the major successes of the 
picture that our universe is one of the many different universes in which low energy physical laws 
take different forms [1] . Such a picture is also suggested theoretically by eternal inflation [2] and the 
string landscape |3]. This elegant picture, however, has been suffering from the predictivity crisis 
caused by an infinite number of events occurring in eternally inflating spacetime. To make physical 
predictions, we need to deal with these infinities and define physically sensible probabilities |1]. 

Recently, a well-defined framework to describe the eternally inflating multiverse has been pro- 
posed based on the principles of quantum mechanics [5]. In this framework, the multiverse is 
described as quantum branching processes viewed from a single "observer" (geodesic), and the 
probabilities are given by a simple Born-like rule applied to the quantum state describing the entire 
multiverse. Any physical questions — either regarding global properties of the universe or outcomes 
of particular experiments — can be answered by using this single probability formula, providing 
complete unification of the eternally inflating multiverse and many worlds in quantum mechanics. 
Moreover, the state describing the multiverse is defined on the "observer's" past light cones bounded 
by (stretched) apparent horizons; namely, consistent description of the entire multiverse is obtained 
in these limited spacetime regions. This leads to a dramatic change of views on spacetime and 
gravity. 

In this paper we present a calculation of the probability distribution of the cosmological con- 
stant in this new framework of the quantum multiverseo We fix other physical parameters and 
ask what values of the cosmological constant A we are likely to observe. In Section [2] we begin 
by reviewing the proposal of Ref. [5], and we then explain how cosmological parameters can be 
calculated in Section El While the framework itself is well-defined, any practical calculation is nec- 
essarily approximate, since we need to model "experimenters" who actually make observations. In 
our context, we need to consider galaxy formation and life evolution therein, which we will do in 
Section |H We present the result of our calculation in Section |5l We find that, in contrast to the 
case with some earlier measures [11] , the measure of Ref. [5] does not lead to unwanted preference 
for a negative cosmological constant — in fact, a positive value is preferred. We find that a simple 
anthropic condition based on metallicity of stars is sufficient to make the calculated distribution 
consistent with the observed value at an order of magnitude level. We conclude in Section [61 

Appendix |X] lists formulae for galaxy formation used in our analysis. Appendix [Bl discusses the 
anthropic condition coming from metallicity of stars. 



^For earlier studies of the cosmological constant in the context of geometric cutoff measures, see Refs. P"[TO]. 



2 Probabilities in the Quantum Mechanical Multiverse 

Here we review aspects of the framework of Ref. [5] which are relevant to our calculation. In 
this framework, the entire multiverse is described as a single quantum state as viewed from a 
single "observer" (geodesic). It allows us to make well-defined predictions in the multiverse (both 
cosmological and terrestrial), based on the principles of quantum mechanics. 

Let us begin by considering a scattering process in usual (non-gravitational) quantum field 
theory. Suppose we collide an electron and a positron, with well-defined momenta and spins: \e^e~) 
at t = — oo. According to the laws of quantum mechanics, the evolution of the state is deterministic. 
In a relativistic regime, however, this evolution does not preserve the particle number or species, so 
we find 

^(t = -oo) = |e+e") -^ ^(t = +00) = Ce I e+e") + c^ I /i+/i~) H , (1) 

when we expand the state in terms of the free theory states (which is appropriate for t — ?■ ±00 when 
interactions are weak). The Hilbert space of the theory is (isomorphic to) the Fock space 

00 
^ = 0^fp, (2) 

71=0 

where "Hip is the single-particle Hilbert space. Various "final states," |e"'"e~) , |yU."'"yU.~) , ■ ■ • , in Eq. ([1]) 
arise simply because the time evolution operator causes "hopping" between different components 
of the Fock space in Eq. ([2]). 

The situation in the multiverse is quite analogous. Suppose the universe was in an eternally 
inflating (quasi-de Sitter) phase S at some early time t = Iq. In general, the evolution of this state 
is not along the axes determined by operators local in spacetime. Therefore, at late times, the state 
is a superposition of different "states" 

^(t = to) = |S) -> *(t) = ^Q(t) l(cosmic) configuration i) , (3) 

i 

when expanded in terms of the states corresponding to definite semi-classical configurations. The 
Hilbert space of the theory is (isomorphic to) 

"H = ^j^'H^^, T-Lm = "Hx^bulk ® "Ha^, horizon, (4) 

M 

where Hm is the Hilbert space for a fixed semi-classical spacetime M., and consists of the bulk 
and horizon parts "HA^.buik and "Ha/j, horizon- (The quantum states are defined on the "observer's" 
past light cones bounded by apparent horizons.) The final state of Eq. ([3]) becomes a superposition 
of different semi-classical configurations because the evolution operator for \l/(t) allows "hopping" 
between different 1-Lm in Eq. (jlj). 



As discussed in detail in Ref. [5], any physical question can be phrased as: "Given what we know 
about our past light cone, A, what is the probability of that light cone to have properties B as 
well?" This probability is given by 

^^^'^^- jdtimioAim) ' ^^ 

assuming that the multiverse is in a pure state |^E'(t)). (The mixed state case can be treated 
similarly.) Here, Oa is the projection operator 

OA = ^\aA,i) {aA,i\y (6) 

i 

where \aA,i) represents a set of orthonormal states in the Hilbert space of Eq. (jl]), i.e. possible past 
light cones, that satisfy condition A (and similarly for Oahb)- Despite the fact that the t integrals 
in Eq. (^ run from t = tg to 00, the resulting P{B\A) is well-defined, since |^I/(t)) is continually 
"diluted" into supersymmetric Minkowski states [3]. 

The formula in Eq. ([5]) (or its mixed state version) can be used to answer questions both regarding 
global properties of the universe and outcomes of particular experiments. This, therefore, provides 
complete unification of the two concepts: the eternally inflating multiverse and many worlds in 
quantum mechanics [5]o To predict/postdict physical parameters x, we need to choose A to select 
the situation for "premeasurement" without conditioning on x. We can then use various different 
values (ranges) of x for B, to obtain the probability distribution P{x). In the next section, we 
discuss this procedure in more detail, in the context of calculating the probability distribution of 
the vacuum energy, x = p\ = A/StiGn. 

3 Predicting/Postdicting Cosmological Parameters 

In order to use Eq. ([5]) to predict/postdict physical parameters, we need to know the relevant 
properties of both the state |^E'(t)) (or its bulk part Pbuik = Tr horizon I'^it)) (^(^)l) ^^^ the operators 
Oa and Oahb- Here we discuss them in turn. 

In general, the state {"^{t)) depends on the dynamics of the multiverse, including the scalar 
potential in the landscape, as well as the initial condition, e.g. at t = to- Given limited current 
theoretical technology, this introduces uncertainties in predicting physical parameters. However, 
there are certain cases in which these uncertainties are under control. Consider x = p\. We are 



^The claim that the multiverse and many worlds are the same has also been made recently in Ref. [12], but the 
physical picture there is very different. Those authors argue that quantum mechanics has operational meaning only 
under the existence of causal horizons because making probabilistic predictions requires decoherence with degrees of 
freedom outside the horizons. Our picture does not require such an extra agent to define probabilities (or quantum 
mechanics). The evolution of our 'l'(i) is deterministic and unitary. 



interested only in a range a few orders of magnitude around pA,obs — (0.0024 eV)^ [13], which is 
tiny compared with the theoretically expected range — Mpj ^ pA ^ Mp^. Therefore, unless the 
multiverse dynamics or initial condition has a special correlation with the value of the vacuum 
energy in the standard model (SM) vacua, we expect that the probabilities of having these vacua in 
|\l/(t)) is statistically uniform in x within the range of interest. (This corresponds to the standard 
assumption of statistical uniformity of the prior distribution of pA [1].) The distribution of x = pA is 
then determined purely by the dynamics inside the SM universes, i.e., the probability of developing 
experimenters who actually make observations of the vacuum energies. 

Let us now turn to the operators Oa and Oahb- In order to predict the value of the vac- 
uum energy which a given experimenter will observe, we need to choose Oa to select a particular 
"premeasurement" situation for that experimenter, i.e. 

„/ X , n/T-,1 ^\ i A '■ a particular "premeasurement" situation ,_. 

P[PA) dpA = P{B\A), < r,. ^ . v^Pnnm PriPrcrv <- n ^ ri n ^ 

I JD . Pa < vacuum energy < pA + "Pa; 

where P{B\A) is defined in Eq. (|5]). Here, we have assumed that the number of SM vacua is suffi- 
ciently large for pA to be treated as continuous in the range of interest. In general, the specification 
of the premeasurement situation can be arbitrarily precise; for example, we can consider a particular 
person taking a particular posture in a particular room, with the tip of the light cone used to define 
|^(t)) located at a particular point in space. In practice, however, we are interested in the vacuum 
energy "a generic observer" will measure. We therefore need to relax the condition we impose as 
A; in other words, we need to "coarse grain" the premeasurement situation. In fact, some coarse 
graining is always necessary when we apply the formalism to postdiction (see discussions in Ref. [5]). 
What condition A should we impose then? To address this issue, let us take the semi-classical 
picture of the framework, discussed in Section 2 of Ref. [5]. In this picture, the probability is given 
by 

P{B\A) = hm ^, (8) 

where N'a is the number of past light cones that satisfy A and are encountered by one of the n 
geodesies emanating from randomly distributed points on the initial hypersurface at t = to- (This is 
equivalent to Eq. (jS]) in the regime where the semi-classical picture is valid.) Since we vary only pA, 
all the SM universes look essentially identical at early times when the vacuum energy is negligible. 
The assumed lack of statistical correlation between pA and the multiverse dynamics then implies 
that we can consider a fixed number of geodesies emanating from a fixed physical volume at an 
early time (e.g. at reheating) in these universes, and see what fraction of these geodesies find the 
"premeasurement" situation A in each of these universes. 

Given that we are focusing on the SM universes in which only the values of pa are different, it 
is reasonable to expect that all the experimenters look essentially identical for different pA, at least 



statistically — in particular, we assume that they have similar sizes, masses, and lifetimes. With this 
"coarse graining," the condition A can be taken, e.g., as: the geodesic intersects with the body of 
an experimenter at some time during their life. In practice, this makes the probability proportional 
to the fraction of a fixed comoving volume at an early time that later intersects with the body of 
an observer. Note that the details of the condition A here do not matter for the final results — for 
example, we can replace the "body" by "head" or "nose" without changing the results because its 
effect drops out from the normalized probability. Thus, in this situation (and any situation in which 
condition A can be formulated entirely in terms of things directly encountered by the geodesic), 
the semi-classical approximation to the scheme of Ref. |5] can be calculated as the fat geodesic 
measure outlined in Ref. [T3]jf| We emphasize that the consistent quantum mechanical solution to 
the measure problem in Ref. [5] forces this choice on us. 

We can now present the formula for -P(pa) in a more manageable form. Since the probability 
for one of the geodesies to intersect an experimenter is proportional to the number of experimenters 
and the density of geodesies, we have 

P(Pa) OC ^ iVobs.a Pgcod,a, (9) 

a 6 habitable galaxies 

where Nohs,a and Pgcod,a are, respectively, the total number of observers/experimenters and the 
density of geodesies in a "habitable" galaxy a. Here, we have approximated that Pgeod.a is constant 
throughout the galaxy a. Note that since we count intersections of experimenters with geodesies, 
rather than just the number of observers (as in much previous work, e.g. [6]), our results differ from 
such previous results by our factor of Pgeod.a- Our remaining task, then, is to come up with a scheme 
that can "model" A''obs,a and Pgcod.a reasonably well so that the final result is not far from the truth. 

4 Approximating Observers 

In this section, we convert Eq. (|9]) into an analytic expression that allows us to compute -P(pa) 
numerically. We focus on presenting the basic logic behind our arguments. Detailed forms of the 
functions appearing below, e.g. F{M, t) and H(t'; M, t), as well as useful fitting functions, are given 
in Appendix \M 

Let us begin with Nobs,a- We assume that, at a given time t, the number of observers arising in 
a given galaxy a is proportional to the total number of baryons in a 






it)OCNB,a{t), (10) 



■^To our knowledge, no detailed study of the probability distribution of the cosmological constant according to the 
fat geodesic measure has been published prior to this work. 



as long as stars are luminous. This assumption is reasonable if the rate of forming observers is 
sufficiently small, which appears to be the case in our universe. To estimate the number of baryons 
existing in galaxies, we use the standard Press-Schechter formalism [15] , which provides the fraction 
of matter collapsed into halos of mass larger than M by time t, F{M, t). Since the amount of baryons 
collapsed is proportional to that of matter, we can use this function F to estimate the number of 



observers and finql 



PiPA) i - jdtjdM ^^^^ Pgcod(M, t). (11) 

The expression of Eq. flTTl) does not take into account the fact that forming intelligent observers 
takes time, or that observers appear only when stars are luminous (which we postulate, motivated 
by the assumption that we are typical observers). To include these effects, we use the extended 
Press-Schechter formalism [16], which gives the probability H{t';M,t) that a halo of mass M at 
time t virialized before t'. The probability density -P(pa) can then be written as 

P{Pa) i - jdt JdM ^^ {Hit - tevoi; M, t) - H{t - tburn; M, t)} Pgeod(M, t), (12) 

where tevoi and tbum are the time needed for intelligent observers to evolve and the characteristic 
lifetime of stars which limits the existence of life, respectively. 

The density of geodesies Pgeod{M, t) is proportional to that of a dark matter halo of mass M at 
time t, which is given by its average virial density: 

P.UM, i) « (4j^ ) / *' ..„«') ^^. (13) 

where t* = min{t, tstop(Af)} with tstop(^) the time after which the number of halos of mass M 
starts decreasing, i.e. when merging into larger structures dominates over formation of new halos: 
d'^F/dMdt\t=tstop{M) = 0. (For the explicit expression of pvir, see Appendix [Aj) In the interest of 
speeding up numerical calculation, we approximate this by the virial density at the time when the 
rate of matter collapsing into a halo of mass M, i.e. —d'^F/dMdt, becomes maximum: 

Pgeod(M,t)f^p,ir(r(M)), (14) 

where t{M) is given by 



d^F{M, t) 



dMdf^ 



= 0. (15) 

t=T{M) 



This approximation is indeed reasonable at the level of precision we are interested in: it works at 
the level of 20% for t ^ 1.7r(M) where the contribution to -P(pa) almost entirely comes from. 



Note that the sign of dF/dM is negative because of the definition of F. 



Finally, there will be several additional anthropic conditions for a halo to be able to host intelli- 
gent observers. For example, the mass of a halo may have to be larger than some critical value M-ann 
to efficiently form stars [17], and smaller than M^ax for the galaxy to be cooled efficiently [HI [T9] . 
Considering these factors, we finally obtain from Eqs. flT2|) and f lT4|) 

P{Pa) = -^r dt dM ] ' ^ {Hit - tevoi; M, t) - H{t - tburn; M, t)} p,ir(r(M)) n{M, t), 

(16) 
where A^ is the normalization factor. Here, 

oo for Pa > 



and we have put anthropic conditions besides Mmjn ^ax in the form of a function n. Note that F, 
H, and pvir (and possibly n) all depend on the value of the vacuum energy p^; see Appendix lAl 

In summary, (dF/dM) {H\t_t^^^^ — H\t-t^^^^^)ndMdt counts the (expected) number of observers 
in halos with mass between M and M + dM at time between t and t + dt, and Pvir(T) is proportional 
to the density of geodesies in such a halo and time, and so Eq. (ITBIl gives the probability by 
counting their intersections (as in Eq. ([9])), with n implementing some anthropic conditions. One 
well-motivated origin for n is metallicity of stars, which affects the rate of planet formation (see 
e.g. Refs. [20l [21]). Here we simply model this effect by multiplying some power m of integrated 
star formation up to time t — tevoi, which we assume to be proportional to the integrated galaxy 
formation rate for M > Mmm: 

n{M, t) oc (F(M„,i„, min{t - tevoi, tstop}) - F{M, min{t - tevoi, tstop}))", (18) 

where tgtop is determined by d{F{M^i^, t') —F{M, t')}/dt'\^,^^^^^ = 0. (For the derivation of this ex- 
pression, see Appendix [HI) Motivated by the observation that the formation rate of certain (though 
not Earth-like) planets is proportional to the second power of host star metallicity [21] , we consider 
the case m = 2, as well as ?7i = ijj 

Another possible anthropic condition comes from the fact that if a halo is too dense, it may not 
host a habitable solar system because of the effects of close encounters [23]. Following Ref. [18], we 
assume this anthropic condition to take the form 

/I 
n^a^v^ < — , (19) 



^ There is observational data for metallicity of galaxies in our universe [32], which our crude model here does not 
reproduce quantitatively. However, when we straightforwardly extrapolate the empirical data to other universes, the 
same regions of the integrand in Eq. (J16p are suppressed/enhanced so that the effect on our calculation is qualitatively 
the same. It must be noted, though, that the strength of the effect may change; e.g. P{pa) for our model with m = 1 
is qualitatively similar to the distribution obtained using the observation-motivated method with m = 3. 



where n^, a^, w^, and tcr are the density of stars, critical "kill" cross section, relative velocity of 
encounters, and some timescale relevant for the condition. Since n* oc pvir, ^'t ~ ^W oc M^/'^pJ-^^ , 
and (j^ and t^r are (expected to be) independent of M and p^jr, this is translated into 

/ / M \2/7\ 

n{M, t) = Q ( p^ax - Pvir(r(M)) [-^^f-) j , (20) 

where Q{x) is the step function (= 1 for x > and = for x < 0), and we have normalized M by 
The value of pmax is highly uncertain. One way to estimate it is to follow Ref. [18] and take 



-3/ Pvir \ 2 T^ir,ym ( ^ V'^f Pvir V^^ 



n^~(lpc) ), at-vrrAu, ^^t ~ ^vir ~ W (v7 — ) ( ) , (21) 

-Pvir,MW^ V '^p ^JWmW^ ^Pvir.MW^ 

where mp is the proton mass, tau — 1.5 x 10® km is the Sun-Earth distance, and Pvir.MW ~ 
2 X 10"^^ g/cm^, Tvir^MW ~ 5 X 10^ K, and Mmw ~ 1 x IO^^Mq are the virial density, virial 
temperature, and mass of the Milky Way galaxy, respectively. Using t^r ~ tevoi = 5 Gyr, Eq. ( 1T^ 
leads to 

p^ax -9X10=^ Pvir,MW (^^) ' ' ~ 3 X 10^22 g/^j^3^ (22) 

V ilamin ^ 

This corresponds to the constraint from direct encounters, i.e. the orbit of a planet being disrupted 
by the passage of a nearby star. There can also be a constraint from indirect encounters: a pass- 
ing star perturbs an Oort cloud in the outer part of the solar system, triggering a lethal comet 
impact [IB]. For a fixed M, this constraint can be about four orders of magnitude stronger than 
Eq. (I22D 

p^ax ~ 3 X 10^26 g/cm^; (23) 

namely, our Milky Way galaxy may lie at the edge of allowed parameter space. 

In our analysis below, we consider either or both of the above conditions Eqs. f lTS]) and fl2Ul) . In 
the real world, there are (almost certainly) more conditions needed for intelligent life to develop. 
However, incorporating these conditions would likely improve the prediction/postdiction for pA- 
In this sense, our analysis may be viewed as a "conservative" assessment for the success of the 
framework, although it is still subject to uncertainties coming from the modeling of observers. 

5 Distribution of the Cosmological Constant 

Our modeling of observers has several parameters which need to be determined phenomenologically: 
Eq. f lTB]l contains Mmin, ^max, ^cvoi, and tbum, while Eq. (I2U|) contains Pmax- We take the "minimum" 
galaxy mass appearing in Eq. (ITB]) to be 

M^i„ = 2 X lO^^Mo, (24) 
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Figure 1: The normalized probability distribution of the vacuum energy P{pa) as a function of pA/PA,ohs- 
The left panel shows P{pa) with the metallicity condition, Eq. (jlSp . with m = (i.e. no condition; dashed, 
blue), m = 1 (dot-dashed, red), and m = 2 (solid, black). The right panel shows P{pa) with the upper 
bound Pinax, Eq. (pO|) . with Pmax = oo (i.e. no constraint; dashed, blue), 6 x 10~^^ g/cm^ (dotted, purple), 
4.5 X 10~^^ g/cm^ (dot-dashed, red), and 3 x 10~^^ g/cm^ (solid, black). 



below which the efficiency of star formation drops abruptly [IT]. For tcvob and tburn? we take them 
approximately to be the age of the Earth and lifetime of the Sun, respectively: 



''evol 



5Gyr, 



^burn 



10 Gyr. 



(25) 



In our analysis below, we use Eqs. (!24|) and fl25|) ; we do not impose the constraint from galaxy 
cooling, i.e. we set Mmax = C)0. While the values of these parameters are highly uncertain, our 
results are not very sensitive to these values. The dependence of our results on them will be 
discussed at the end of this section. 

In Fig. [H we present the normalized probability distribution for the vacuum energy -P(pa) as a 
function of pA/pA,ohs, under several assumptions about the function n: 

(i) "minimal" anthropic condition: n{M, t) = 1 

(ii) metallicity condition: Eq. (fT8|) with 771 = 1 and 2 

(iii) maximum virial density condition: Eq. (120|) with p^ax = {3xl0~^^, 4.5x10"^^, 6x10"^^} g/cm^, 
which are {1, 1.5, 2} times the value in Eq. fj23|l . 

(The result with Pmax given by Eq. ( 122|) is virtually identical to the case with the minimal anthropic 
condition.) The left panel presents the effects of metallicity, showing (i) and (ii), while the right 
panel those of Pmaxj with (i) and (iii). 

Interestingly, in all cases, our predictions prefer a positive cosmological constant over a negative 
one, as opposed to the situation in earlier measure proposals where strong preferences to negative 
values have been found [TT]. In Tabled] we provide the probabilities of having pA > (and < 0) in all 
six anthropic scenarios. The absence of an unwanted preference towards negative pa is satisfactory. 





P(PA > 0) 


P{PK < 0) 


No condition 


97% 


3% 


Metallicity, m = 1 


87% 


13% 


Metallicity, m = 2 


75% 


25% 


Pmax = 6 X 10-26 g/cm^ 


92% 


8% 


Pmax = 4.5 X 10-26 g/cm^ 


83% 


17% 


Pmax = 3 X 10-26 g/cm^ 


63% 


37% 



Table 1: The probability of observing a positive and negative cosmological constant, P[p\ > 0) and 
P{p!v < 0), for six assumptions on the anthropic condition. In all cases, a positive value is preferred over 
a negative one, consistent with observation. 
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Figure 2: Same as Fig. [H but the horizontal axis now in logarithmic scale. To show the probability density 
per tenfold, the vertical axis is chosen to be pkP{pa)/ P\,ohs- The distributions are normalized in the region 
PA>0. 

especially given that the measure of Ref. |5] was not devised to cure this problem. It comes from 
the fact that the present measure does not have a large volume effect associated with the global 
geometry of anti-de Sitter space, which was responsible for a strong preference for negative pa in 
earlier, geometric cutoff measures [H]. In contrast with these measures, the quantum measure 
of Ref. [S] does not count the number of events; rather, it gives quantum mechanical weights for 
"situations," i.e. quantum mechanical states as described from the viewpoint of a single observer 
(geodesic). The preference towards a positive value comes from the fact that for pA > some 
observers still form after vacuum energy domination, while for pA < it is not possible due to the 
big crunch. 

Figure [1] shows that -P(pa) is always peaked near pA = 0, with the distribution becoming wider 
as the anthropic condition gets weaker. In Fig. [21 we plot the same distributions in logarithmic 
scale for pA/PA,obs, limiting ourselves to pa > 0. To show the probability density per tenfold, the 
vertical axis is chosen as PA-P(PA)/PA,obs- From these figures, we find that our anthropic assumptions 
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Figure 3: The normalized probability distribution P{pk) with the metallicity condition, Eq. (jl8|) . with 
771 = 2. In the upper-left panel, Mmin is varied as 2 x IQ^^Mq (solid, black), 6 x IO^^Mq (dot-dashed, red), 
and 0.67 x 1O^^M0 (dashed, blue); and in the upper-right, Mmax as oo (solid, black), lO^^M© (dashed, 
blue), IQi^^Mq (dot-dashed, red), and 2 x 10^"^ Mq (dotted, purple). The lower left and right panels vary 
tcvoi and tburn as {(solid, black), (dot-dashed, red), (dashed, blue)} = {5,1,0} Gyr and {10,7,15} Gyr, 
respectively. 



lead to results that are consistent with the observed value within one or two orders of magnitude. 
In particular, metallicity alone is enough to bring the agreement to an order of magnitude level. 
This is because mergers, which lead to an increase in metallicity, are suppressed for larger values of 
Pa due to earlier vacuum energy domination. This result is comfortable, especially given that the 
constraint from encounters is effective only if Pmax is close to the Milky Way value, as in Eq. 
Given our crude treatment of observers, we consider these results quite successful. 



Finally, we discuss the sensitivity of our results to variations of M^ 



M„ 



tevoi; and tbu 



which can be thought of as "systematic effects" of our analysis. In Fig. [3l we show the distributions 
of -P(pa) with the m = 2 metallicity constraint, varying the values of Mmm, Mmax, ^evoi, and tburn, 
respectively. We find that, while the detailed shape of -P(pa) does change, our main conclusions are 
robust: (i) There is no strong preference to a negative vacuum energy; in fact, a positive value is 
preferred, (ii) The predicted distribution of pa is consistent with the observed value at an order of 
magnitude level with the metallicity constraint. 
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Figure 4: The normalized probability distribution P{pa) with a metallicity condition: Eq. (llSp with m = 2. 
The light and dark shaded regions indicate those between 1 and 2a, and outside 2a, respectively. The 
observed value /OA/pA,obs = 1 (denoted by a vertical line) is consistent with the distribution at the la level. 

6 Conclusions 

In this paper, we have studied the probability distribution of the cosinological constant (or the 
vacuum energy pa) in the multiverse, using the quantum measure proposed in Ref. [5]. We have 
found that this measure does not lead to a strong preference for negative p\, as opposed to earlier 
measures proposed based on geometric cutoffs, because it does not experience a large volume effect 
associated with the global geometry of anti-de Sitter space. Moreover, we have found that a positive 
value of Pa is preferred, consistent with observation. 

We have found that a simple, intuitive condition based on metallicity is enough to reproduce the 
observed value of p^ at an order of magnitude level. This is comfortable because effects from other 
possible anthropic conditions, such as the ones from encounters, are much more sensitive to the 
details of the conditions. In Fig. |U we present the normalized distribution -P(pa) with the m = 2 
metallicity constraint, where the 1 and 2a regions are indicated. We find that the observed value is 
consistent with the calculated distribution at the la level. 

It would be interesting to refine our analysis including more detailed anthropic effects, such as 
those of star formation. Another possible extension of the analysis is to vary other cosmological 
parameters, such as the primordial density contrast Q and spatial curvature Q^ (at a specified time), 
in addition to p\. We plan to study these issues in the future. 
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V2a{M,t)J ' 



A Press-Schechter Formalism and Fitting Functions 

The Press-Schechter function F is 

F(M, t) = erfc 

where 6c(t) and a{M,t) are given by [TH |2] 
with tcrunch defined in Eq. flT7|) . and 



1.629 + 0.057 e-2-3GivPAi2 
1.686 + 0.165 (-^' 

\ ^crunch 



+ 0.149 
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for Pa > 
for Pa < 0, 



a{M,t) -Qs{M)G{t). 
Here, Q is the primordial density contrast, 

s{M) ~ [(9.1 p-2/^)-°-2^ + {50.5 logio(834 + p-^/^) - 92}' 
where p = M^'^G^ with ^ = Pmattcr/'^7 — 3.7 eV, and 



■0.27 



-1/0.27 



0.206^ 



Pa 



G{t) - < 



tanh2/^(|i7At){l - tanh^-2^(|/7At)} 
+ 1.437(1 - cosh-^/^(|i7At)} 



0.82 



0.549 e'/^G^v^'t^/s 



1 + 0.37 



t 



''Crunch 



2.18 



t 



^crunch 



for Pa > 

for Pa < 0, 



where Ha = ^/SttGnIpaIJ^. 

The function H in the extended Press-Schechter formalism is given by [T6l [25] 



where 



/3(Mi,ti,M2,t2) =erfc 



Uv2(s(Mi)2-s(M2)2) 



VG(ti) G{t2)J 



(26) 



(27) 



(28) 



(29) 



(30) 



(31) 



(32) 
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with s(M) and G{t) defined in Eqs. (129|) and (|30|). 

The virial density as a function of time can be fit, following Refs. [T8|l25]. as the density evolution 
of a closed universe, according to Birkhoff' s theorem. The virial density is then given in terms of 
the density at turn-around rescaled by the ratio of the volumes, pvir = (-Rvir/-Rturn) Upturn- Here, 
Rvir/Rtnrn — ?" 2 at early times (t ^ ^/H\) as well as for |pa| — j- at any fixed t. For positive 
Pa, -RvirZ-Rturn = 2/(-\/3 — 1) ~ 2.73 at late times |26], while for negative pA, Rvn/ Rtum — ^ 2^/^ for 
t -^ ^crunch- Our fit is given by 



P^Ut)^{ r ..J^l ^ ^ (33) 

187rVo.atte.(t)^^;^|^) ':^f_^ . for PA < 0, 

1,2-nAi; / 123.6+7(e ■ *crunch_lj 

where Pmatter is the matter energy density. This fit is accurate to better than ^ 5% and 2% for 
Pa > and < 0, respectivelyjj 

Finally, the time at which most of galaxies of mass M forms, i.e. the solution to Eq. (ITSll . is well 
approximated by the following fitting function: 

3/2 ( 1.880 + ci(a) M + cs{a) M'^ + c^ia) M^ for - 10 ^ ^f^ < 
r(M)/Gyr~%| J c[,(a) + c'i(a)M + 4(a)M3 + c'5(a)M5 for < -^^ <'lO (34) 



g3/2 



c[,'(a) + c'lia) M + 4(a) M^ + 4(a) M^ for 10 < ^^ < 100 



PA,obs 



where a = (pA/PA,obs)(<5obs/Q) , M = log^o 2xio"Mq ' ^^^ 

ci(x) = -0.311 + 1.276 e°-*^27x ^ ^4^2 logroll + |x|°-^}, 

csjx) = 0.470 -0.656e°-^s^- 0.317 logio(0.2 + |x|), 

cgja;) = -0.0142 + 0.0381 e°-^^' + 0.00822 logio(0.05+|a;|), 

4{x) = 1.880- 0.00205 X, 

4{x) = 0.408 + 0.569e-i-oi^' + 0.295 logio(l + x), 

4(x) = 0.277- 0.251 e"^- 0.125 logio(l + a:), (35) 

4lx) = -0.000889 + 0.0151 e"^- 0.00220 logio(l+a;), 

4{x) = 1.880- 0.00205 X, 



4{x) = 0.767 -0.00293 a; - 230 a;-^ 

4lx) = -0.530 + 0.000336 a: + 0.847 x-°-\ 

4{x) = 0.106 -0.0000118a;- 0.125 a;-°-i- 0.0131 logio(-5 + x). 



^For PA < 0, the approximation leading to Eq. (pS)) . i.e. pvir ^ a few (pmattcr + PA), breaks down for t/tcrunch ^ 0.8, 
where we should rather use pvir = (since there is no stable structure forming). However, since pvir in Eq. ([33]) is 
small there anyway, using it up to t/tcrunch = 1 does not lead to a significant error. 
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This fit is accurate to better than ^ 5% for M ^ IO^^Mq (but it becomes worse for smaller M, e.g., 
the accuracy is ~ 12% at M ~ 6 x lO^^M©). For 100 < pA/pA,obs < 150, we use the last expression 
of Eq. (l34l) . good to the level of ~ 10%; and for pA/pA,ohs > 150, we use 

Pvir(r(M)) /(lO-^^ g/cm^) ^ (3.66 + 0.032 a) - (1.36 + 0.0013 a)M, (36) 

which is accurate to the level of ~ 10% up to pA/pA,obs ~ 4500. 

B Anthropic Condition from Metallicity 

In this appendix, we derive the function n arising from the metallicity constraint, Eq. (fTSjl . Suppose 
that in a merging tree of a galaxy j at time t, j is found to have progenitor galaxies i = 1, 2, ■ ■ • 
with varying masses Mi at time t' < t. Note that this also includes accretion, i.e. matter that was 
not part of galaxies of appreciable size, since F{M = 0,t) = 1 in the Press-Schechter formalism, 
where accretion is treated as mergers of extremely tiny galaxies with a large galaxy. 

Now, let us assume that the relative mass fraction in j that came from i and i' is given by 
— dM- I — (Im' ' ^■^' ^^^ ratio of total amount of baryons at time t' in galaxies of type i and i\ 
respectively. This is true within the Press-Schechter formalism as long as Mj j/ -C M^, since then 
the overdensities within spherical top-hat regions containing masses Mj j' and Mj are independent 
of each other at early times. Once Mj^j/ ^ Mj, the assumption is not justified, but in these regimes, 
there can only be a small amount of merging occurring from i,i' to j, implying little contribution 
to metallicity. The assumption, therefore, provides a good approximation. 

Let Xi {i = 1,2, ■ ■ ■) be the fraction of baryons in the universe that formed stars in halos of 
mass Mi at time t'. In our simple model, the star formation rate is proportional to the rate of halo 
formation for masses M > M^i^ and otherwise zero: dxi/dt' oc 0(Mj — Mmin) cPF{Mi,t') / dMidt' . 
The increase in total metal content summed over galaxies of mass Mj is taken to be proportional 
to the star formation rate therein, dxi/dt, so the increase in (linear) metallicity dZi is 

,z.(t') oc ('-^i^Y'^M' « ('-E(M^Yf3Miine(M, - m„,„)< (37) 

^ ' \ dMi J dt' V dMi J dMidV ^ ' ^ ' 

i.e. the total increase in metal content divided by the total mass. The increase in metallicity of 
galaxy j due to stars at time t', then, has to be weighted by the relative matter fraction of galaxies 
i, as described above: 



i^i dM, ^^^y^ ) 2^i dNhdt' ^y^^^^ ^^^^ 

— rx" 

dF{Mi,t') Y^ dF{Mi,t') 



A7 (+l\ — ^^» dMj "'^ty^ / L^i dMjdt' ^V'-'^ "--'mm; , 

""^jV')— ^ dF{Mi,t') ^ y^ dF(Mi,t ') ^'^ ' ^'^^ 

^i dMi ^i dMi 



where we must normalize to the total mass of galaxy j at each time t' . 
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In the continuum limit, the sum over i becomes an integral over masses. Therefore, the metal- 
licity of galaxy j of mass M at time t is 

, , /•* , L^'dM' '"'J/i!';/'^ e(M' - M^in) /■* , MF(M^in,t') - F(M,t')} , , 







J^^^dM''-^^ Jo l-F{M,t' 



where t = min{t, tstop}- The constraint that one cannot accumulate negative metallicity determines 
the timescale tstop as a solution to 



^{F{M^,^,t') - F{M,t')} 



0, (40) 



t'=t 



stop 



at which time merging of galaxies of mass Mmin < M' < M into those more massive than M begins to 
dominate over formation of new galaxies in this mass region. Since merging into larger structures is 
also occurring at earlier times, one expects that we slightly underestimate the metallicity. However, 
in practice, the formation of new galaxies in this mass range and mergers into structures beyond 
are well separated in time, so a simple cutoff at tstop is sufficient. 

Now, since F{M,t') < F(M,tstop) = erfc(l/A/2) ~ 0.317, the denominator of Eq. §9^ is always 
between 0.68 and 1; in fact, it is very close to 1 in most of the parameter regions. Therefore, we 
can safely ignore the denominator of Eq. fl39|) and obtain 

ZiM, t) cc {F{M^^, i) - F{M, i)) . (41) 

In general, the probability of forming planets is expected to be proportional to some power m of 
the metallicity [201 EI] • This gives 

n(M,t) = Z(M,t-tevoi)'", (42) 

which is Eq. (ITSl) in the text. 



References 

[1] S. Weinberg, Phys. Rev. Lett. 59, 2607 (1987). 

[2] A. H. Guth and E. J. Weinberg, Nucl. Phys. B 212, 321 (1983); A. Vilenkin, Phys. Rev. D 27, 
2848 (1983); A. D. Linde, Phys. Lett. B 175, 395 (1986); Mod. Phys. Lett. A 1, 81 (1986). 

[3] R. Bousso and J. Polchinski, JHEP 06, 006 (2000) |arXiv:hep-th/0004134| ; S. Kachru, 
R. Kallosh, A. Linde and S. P. Trivedi, Phys. Rev. D 68, 046005 (2003) 
|arXiv:hep-th/0301240| ; L. Susskind, |arXiv:hep-th/0302219[ M. R. Douglas, JHEP 05, 046 



(2003) [arXiv:hep-th/0303194| . 



16 



[4] For reviews, see e.g. A. H. Guth, Phys. Kept. 333, 555 (2000) |arXiv:astro-ph/0002156] ; 
A. Vilenkin, J. Phys. A 40, 6777 (2007) |arXiv:hep-th/0609193| ; S. Winitzki, Lect. Notes 
Phys. 738, 157 (2008) |arXiv:gr-qc/0612164| ; A. Linde, Lect. Notes Phys. 738, 1 (2008) 
jarXiv:0705.0l64l [hep-th]]. 

Y. Nomura. larXiv: 1104. 2321 [hep-th]. 

H. Martel, P. R. Shapiro and S. Weinberg, Astrophys. J. 492, 29 (1998) 
|arXiv:astro-ph/9701099| ; G. Efstathiou, Mon. Not. Roy. Astron. Soc. 274, L73 (1995). 

J. Garriga, M. Livio and A. Vilenkin, Phys. Rev. D 61, 023503 (2000) |arXiv: astro-ph/ 99062 1 0] ; 
L. Pogosian and A. Vilenkin, JCAP 01, 025 (2007) |arXiv:astro-ph/0611573| . 

R. Bousso, R. Harnik, G. D. Kribs and G. Perez, Phys. Rev. D 76, 043513 (2007) 
|arXiv:hep-th70702115j ; J. M. Cline, A. R. Frey and G. Holder, Phys. Rev. D 77, 063520 
(2008) |arXiv:0709.4443] [hep-th]]. 

A. De Simone, A. H. Guth, M. P. Salem and A. Vilenkin, Phys. Rev. D 78, 063520 (2008) 



[5] 
[6] 

[7 

[8 

[9 
[10 

[11 

[12 
[13 

[14 

[15 
[16 
[17 

[18 

[19 



|arXiv:0805.2173i [hep-th]]. 

R. Bousso, L. J. Hall and Y. Nomura, Phys. Rev. D 80, 063510 (2009) jarXiv: 0902. 2263) [hep- 
th]]. 

M. P. Salem, Phys. Rev. D 80, 023502 (2009) |arXiv:0902.4485] [hep-th]]: R. Bousso and S. Le- 
ichenauer, Phys. Rev. D 81, 063524 (2010) |arXiv:0907.4"9T7l [hep-th]] : R. Bousso, B. Freivogel, 
S. Leichenauer and V. Rosenhaus. larXiv: 1012.28691 [hep-th]. 

R. Bousso and L. Susskind, larXiv: 1105.3796 1 [hep-th]. 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) [arXiv: 1001. 45381 
[astro-ph.CO]]. 

R. Bousso, B. Freivogel and I. S. Yang, Phys. Rev. D 79, 063513 (2009) [ arXiv:0808.3770l 

[hep-th]]. 

W. H. Press and P. Schechter, Astrophys. J. 187, 425 (1974). 

C. Lacey and S. Cole, Mon. Not. Roy. Astron. Soc. 262, 627 (1993). 

G. Kauffmann et al. [SDSS Collaboration], Mon. Not. Roy. Astron. Soc. 341, 54 (2003) 
|arXiv:astro-ph/0205070| . 

M. Tegmark, A. Aguirre, M. J. Rees and F. Wilczek, Phys. Rev. D 73, 023505 (2006) 
[arXiv: astro-ph/05 1 1 774] . 

A. Cattaneo, A. Dekel, J. Devriendt, B. Guiderdoni and J. Blaizot, Mon. Not. Roy. Astron. 
Soc. 370, 1651 (2006) [arXiv:astro-ph/0601295| . 



17 



[20] C. H. Lineweaver, Icarus 151, 307 (2001) |arXiv:astro-ph/0012399] . 

[21] D. A. Fischer and J. Valenti, Astrophys. J. 622, 1102 (2005); D. Grether and C. H. Lineweaver, 
Astrophys. J. 669, 1220 (2007). 

[22] R. Maiolino et al, Astron. Astrophys. 488, 463 (2008) |arXiv:0806.24T0l [astro-ph]]; 
P. S. Behroozi, C. Conroy and R. H. Wechsler, Astrophys. J. 717, 379 (2010) [arXiv:1001.0015 

[astro-ph. CO]]. 

[23] M. Tegmark and M. J. Rees, Astrophys. J. 499, 526 (1998) |arXiv:astro-ph/9709058| . 

[24] J. A. Peacock, Mon. Not. Roy. Astron. Soc. 379, 1067 (2007) jarXiv:0705.0898l [astro-ph]]. 

[25] R. Bousso and S. Leichenauer, Phys. Rev. D 79, 063506 (2009) |arXiv:0810.3"044l [astro-ph]]. 

[26] O. Lahav, P. B. Lilje, J. R. Primack and M. J. Rees, Mon. Not. Roy. Astron. Soc. 251, 128 
(1991). 



